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The  problem  of  receiving  intelligence  through  a medium  whose 
properties  are  time  varying  has  become  an  important  one  recently  as 
communication  techniques  are  extended  to  the  upper  atmosphere,  or  trop- 
o sphere,  and  to  the  underwater  domain  of  submarines.  Both  the  tropo- 
sphere and  sea  water  are  time  varying  in  a number  of  ways,  and  the  time 
variations  have  a wide  range  of  characteristic  periods.  Some  of  these  are 
longs  such  as  die  diurnal  and  seasonal  variations  of  underwater  sound  ve- 
locity profiles.  Others  are  quite  short  and  may  result  from 
storms  and  solar  flares„  which  affect  the  ion  density  of  ths  troposphere, 
or  the  wave  motion  of  the  seawater  surface  and  the  relative  motion  of 
different  underlying  layers  of  water.  This  paper  JS^B  be  concerned  with 
tuch  short-time  temporal  instabilities,  namely,  those  that  ha  >e  an  oppor 
tunity  to  distort  the  time  base  of  a transmitted  or  received  signal  (wave- 
form} during  the  time  duration  of  the  coded  signal.  Such  distortions,  even 
though  slight,  can  be  expected  to  play  havoc  with  the  usual  type  of  matched 
filter  or  correlator  receiver*  which  is  actually  an  optimum  receiver  when 
temporal  distortions  are  not  present^  Moreover,  as  signal  codes  of 
longer  duration  are  used  in  an  attempt  to  put  more  energy  into  the  medium 
and  attain  more  processing  gain  in  the  receiver,  the  liklihood  of  temporal 
distortions  occurring  within  the  time  duration  of  the  signal  code  becomes 
more  pronounced.  It  will  be  shown  that  the  optimum  receiver  in  the  pres- 
ence of  temporal  distortions  is  a particular  form  of  an  adaptive  receiver. 

The  communication  problem  treated  here  is  considered  in  the 
general  sense;  it  includes  the  transmission  of  a predetermined  waveform 
coded  in  any  manner,  and  the  reception  of  a temporally  distorted  version 
of  this  waveform  plus  additive  white  Gaussian  noise.  Thus,  it  pertains  to 
message  transmissions  for  purposes  of  sending  intelligence  from  one  point 
to  another  as  well  as  to  radar  and  sonar  systems  designed  principally  for 
detection.  Actually  the  detection  aspect  of  an  optimally  encoded  waveform 


* P,  M.  Woodward,  Probability  and  Information  Theory,  with  Applica- 
tions to  Radar  (Pci-gsmon  Fbress,  Oxford,  Chap,  4. 
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will  be  emphasized  in  this  paper;  no  attempt  will  be  made  to  treat  the 
choice  of  codes  or  the  best  utilization  of  channel  capacities  in  a Shannon 
sense, 

A second  important  cause  of  temporal  incoherence,  in  addition 
to  the  vagaries  of  the  medium,  becomes  apparent  when  we  consider  target 
detection  by  reflection  of  a long  duration  encoded  signal  A target,  to  be 
worthy  of  its  mission,  obviously  has  the  capability  of  changing  its  location 
m time  and  of  accelerating,  and  most  targets  can  do  both  during  the  time 
they  are  illuminated  by  a radar  or  sonar  transmitter.  The  effect  of  such 
target  motion  is  to  multiply  the  time  base  of  the  reflected  waveform  by  a 
constant  scale  factor,  in  the  case  of  uniform  relative  motion,  or  by  a 
variable  scale  factor,  in  the  case  of  accelerated  motion.  The  former 
effect  produces  the  Doppler-distorted  return  signal  for  which  satisfactory 
techniques  have  been  developed.  The  latter  effect  appears  as  an  unwanted 
temporal  distortion  of  the  signal  a '’rubberizing"  of  the  time  base  which  is 
quite  analogous  to  that  imposed  by  the  time  variations,,  e,  g. , ocean  surface 
motion,,  of  the  medium. 

An  assumption  will  be  made  that  signals  arriving  at  the  receiver 
via  several  different  paths  in  the  medium  are  resolvable  even  though  over- 
lapping, In  other  words,  the  various  multipath  arrivals  are  separated  in 
time  so  that  the  output  of  a correlation  receiver  is  a set  of  resolved  corre- 
lation peaks  corresponding  to  the  various  multipaths.  Part  of  a proper  en- 
coding procedure  is  the  choice  of  waveform  which  possesses  a good  ambigu- 
ity function  and  hence  good  time  and  Doppler  resolution,  ^ and  therefore  good 
multipath  resolution.  Moreover,  the  physical  nature  of  the  medium  is  often 
such,  particularly  in  sound  propagation  through  sea  water,  that  the  various 
arrivals  are  separated  in  time  by  appreciable  fractions  of  a second,  more 
than  enough  to  resolve  an  optimally  coded  waveform  of  wide  bandwidth. 

The  well  known  optimum  detector  for  a fully  coherent  signal  cor- 
rupted by  additive  Gaussian  noise  is  one  which  cross  correlates  the  received 
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signal  with  a pre  recorded  version  of  th©  transmitted  signal  This  result 
can  be  derived  from  a likiihood  function  as  Woodward  and  others  have 
shown  As  a next  step  to  increasing  the  maximum  range  of  detection,  the 
transmitted  waveform  is  lengthened  (in  the  case  of  a peak-power  limited 
transmitter*  no  as  to  couple  more  energy  into  the  medium  and  to  increase 
the  frequency*  or  Doppler  velocity,  resolution  for  a moving  target.  How 
ever,  in  the  presence  of  the  temporal  instabilities  already  discussed*  ex- 
tending the  time  duration  of  the  signal  becomes  self-  defeating  after  a point 
Despite  the  increase  in  input  signal  to  noise  ratio  obtained  thereby,  the 
greater  the*  frequency  resolution  (resolution  spacing  - 1/T»  possessed  by  a 
waveform  by  reason  of  t&  time  duration  T the  more  susceptible  is  its 
autocorrelation  function  to  slight  distortions  of  the  time  base  Thus*  it  is 
quite  possible  to  "overdesign' : a detection  system  by  inadvertently  building 
so  much  resolution  into  it  that  is  tails  to  recognize  a slightly  distorted,  tar 
get  reflected  signal  in  spite  of  an  adequately  high  input  signal  to  noise  ratio 
fa  fact,  such  an  "overdesigned"  system  can  easily  fail  to  detect  signals 
which  would  have  been  recognized  unambiguously  by  simpler,  less  costly 
receivers  which  processed  the  received,  energy  coherently  over  a shorter 
time  span  Thus  is  ^nvu>us  that  an  optimum  detector  should  be  matched 
somehow  to  the  "coherence  times"  of  the  medium  and  any  moving  targets 
which  may  be  present.  It  is  less  obvious  exactly  what  form  an  optimum  re 

reiver*  so  matched*  should  take,  although  this  problem  has  been  treated 

. 4 5 

in  the  radar  case  for  a received  signal  distorted  by  multiplicative  noise.  " 

For  certain  sonar  applications  particularly,  it  is  more  pleasing  and  specific 

to  treat  the  exact  case  of  a time  base  distortion.  The  following  paragraphs 

derive  the  form  of  an  optimum  receiver  from  the  standpoint  of  a maximum 

likiihood  detector  when  time  base  distortions  are  present. 

We  will  consider  a received  signal  f(t  which  is  comprised  of 

additive  white  Gaussian  noise  n(t«  and  a distorted  version  sfk  t i - T(t.l 

1 o o J 
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C.  W.  Hel  stroro,,  Statistical  Theory  of  Signal  Detection  (Per  gam  on  Press, 
Oxford,  I960!,  p,  64,  73 

4 R,  Rojas,,  private  communication, 

* P.  Bello,  "Joint  estimation  of  delay,  Doppler*  and  Doppler  rate,  M IRE 
Trans, , P„G„IT0,  IT-6,  330  (1960). 
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of  the  transmitted  signal  s(t)  r 


f(ti  niti  + r8[kQt  » - T(tj]  , 


(II 


The  delay  parameter  t represents  the  round  trip  time  delay  for  propaga- 
tion out  to  the  target  and  back,  averaged  over  the  signal  duration  The 
term  T|t>  introduces  the  time  base  distortion,,  which  in  itself  is  a function 
of  time  The  term.,  a constant  accounts  for  the  time  base  compression 
or  expansion  resulting  from  the  target  s average  speed  relative  to  the 
transmitter- receiver  location.  The  quantity  r is  a real  constant  which 
represents  attenuation  of  the  transmitted  signal  The  r might  more  gen 
e rally  have  been  considered  complex  and  time  varying,  in  which  case  its 
phase  would  represent  a phase  shift  of  the  signal  carrier  and  the  time  flue- 
tuation  < ou!d  account  for  scintillation  effects.  While  such  considerations 
are  applicable  to  radar„  it  is  the  contention  here  that  phase  shifts  of  the 
carrier  do  not  adequately  represent  the  usual  time  base  distortions  in  un- 
der* ater  ». gnat  propagation 


Only  statistical  knowledge  is  available  to  the  observer  for  n(ti  and 
T(t  in  Eq  ( 1 Their  forms  as  explicit  functions  Of  t are  not  known,,  but 
they  can  be  described  in  terms  of  the  probability  distributions  of  their  ampli- 
tudes 9 


(2> 


and 
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for  any  time  t „ 
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The  probability  distribution  upon  *f  is  assumed  to  be  Gaussian 
not  only  because  of  its  simplicity  and  its  frequent  close  approximation  to 
actual  probability  distributions  in  nature,  but  because  by  the  central  limit 
theorem  the  actual  distribution  will  tend  to  a Gausrian  shape  if  a number 
of  independent  and  random  factors  contribute  to  the  r(t)  distortion  at 
any  given  time.  In  tropospheric  and  underwater  propagation  over  long 
distances  this  is  the  case,  for  the  resultant  time  base  distortion  is  a cas- 
cading of  many  spatially  separated  distorting  effects,  For  example,  a 
principal  mode  for  transmission  of  underwater  sound  over  long  distances 
is  the  RSK  mode  (refraction  surface  reflection)  which  entails  multiple 
surface  or  near  surface  reflections  before  arrival  at  the  receiver.  The 
sea  wave  amplitudes  and  phases  at  the  spatially  separated  surface  reflec- 
tions can  be  considered  independent,  and  reflection  from  each  of  the 
moving  surfaces  acts  sequentially  as  a time  base  distorter  upon  the  prop- 
agated signal  Thus,  the  resultant  distortion  can  be  thought  of  as  an  op- 
erator which  is  the  product  of  a number  of  other  independent  operator  s 
cascaded  to  produce  the  resultant  effect, 

£quations  (2s  and  (3)  g:.ve  the  probability  distributions  at  any  time 
t „ but  the  joint  probability  distributions  for  all  times  included  in  the  sig- 
nal duration  are  actually  of  interest,  For  this  it  is  convenient  to  use 
Shannon  s concept  of  number  of  'degrees  of  freedom"  in  a waveform. 
Shannon  showed  that  a waveform  of  duration  T containing  frequencies  with 
a bandwidth  W has  2WT  degrees  of  freedom.  In  other  words,  any  such 
waveform  can  be  fully  specified  by  a set  of  real  time  samples  j W sec  apart 
taken  over  the  time  span  T , or  by  a set  of  complex  frequency  components 
spaced  IT  cps  apart  summed  over  the  bandwidth  W,  Since  the  samples 
are  statistically  independent,  the  joint  probability  distribution  for  any 
arbitrary  noise  sample  can  be  written  as  a multivariate  normal  distribu  - 
tion without  cross  correlation  terms: 


Wq  [n(t>]  * 
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Equation  (5  ran  be  rewritten  more  compactly  by  considering  the  2WT 
independent  components  t to  form  a vector  in  a signal  space  of  2WT 
orthogonal  dimensions  Then  Kc  {5]  becomes 
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v re  the  arrow  above  a quantity  indicates  a vec  tor.  For  convenience, 
we  wril«  describe  t Eq„  (3  „ in.  the  same  space  of  AWT  dimeneiors 

even  though  the  number  of  independent  components  in  the  vec  tor  t is 
ZWJT  _ where  is  the  bandwidth  of  frequencies  involved  in  ■jjti  „ 

Since  >U  arises  in  tne  first  place  because  of  physical  motion  of  the  prop* 
a gating  medium  and  the  target,,  its  frequencies  must  he  band  limited  to  a 
low  pass  region  descriptive  of  these  phenomena  Hence,  it  is  reasonable 
to  assume  that  Wt  « W the  signal  bandwidth...  Then  a description  of 
r in  a space  of  2WT  dimensions  will  require  cross  correlation  coefficients 


2WT 


where  the  diagonal  elements  are  all  one„  and,  as  shown  by  the  Gaussian 
shaped  curve  at  the  center  oi  the  matrix,  the  slightly  off-diagonal  ele- 
ments will  be  less  than  umty„  falling  to  zero  as  their  separation  from 
the  main  diagonal  increases,  This  form  of  the  correlation  matrix  shows 
that  closely  spaced  sample  points  t are  highly  correlated*  with  the 
correlation  dropping  to  zero  as  the  tune  between  sample  points  increases. 
The  shape  of  this  decay  depends  upon  the  time  constants  of  the  physical 
causes  of  the  temporal  distortions;  in  the  case  of  a sonar  system*  it 

would  depend  upon  the  sea  state  and  acceleration  capabilities  of  the  targets 

2 

Finally.,  the  variant  e factor  c ^ is  unchanged  whether  Wq(t,  is  described 
in  a space  of  2WMT  or  2WT  dimensions.  The  variance  is  indicative  of  the 
range  of  amplitude  excursions  in  Tftl  and  also  depends  upon  sea  conditions 

target  acceleration  capability,  but  it  it  independent  of  the  shape  of  die 
frequency  spectrum  «£  tftl.  Rewriting  Eq»  O)  as  a multivariate  normal  dis- 
tribution in  a apace  of  2WT  dimensions,,  we  take  account  of  the  cross 
correlation  of  the  t.  terms  as  follows; 

r (pict 
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2cr  r lp  I 


where 


ipiC  • cofactor  of  the  correlation  matrix,  and 
| p ■ determinant  of  the  correlation,  matrix 
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To  iefine  «n  optimum  receiver  #e  Woodward' a conceptf 

of  a ma  tjoium  (ikli!  ood  detector  * Certainly  the  beat  that  any  receiver 

could  provids  w ould  be  the  a posteriori  probability  distribution 

Wt'r  V T I?  ? of  T * and  V „ git  en  th<£  recel  1 aifflca!  T „ whese 
o o o w 

W|  *•  V rolf:  ‘ WOl?'  V V L ‘T  ]T-  kO  To' 
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No  further  information  could  be  obtained  or  desired  from  the  received 
signal  f 

In  Eq  (&!  c is  a constant  and  W (r  , k , t is  the  a priori  proh 
ability  densitv  di etr  buticns  of  T»k  and  r which  can  be  factored  ae - 
cording  to  the  lower  eorpression  in  Eq  (3  into  two  independent  terms 
W fl  !■  o a function  of  target  motion  (position  and  relative  speedj » a-ud 
W h ) , which  includes  the  perturbing  effects  oi  the  medium  and  accelera 
tl«m  effects  She  -a  tss  tne  time  T • L.|?  |f  ,,k  . t } is  the  iikli 

h"  . ::u.  son  which  is  proportional  to  the  conditional  probability 
Wlf  ir  „ k g t > cons  de red  as  a function  of  rr  s i.  and  r The  objective 
oi  ar  optimum  receiver  is  to  evaluate  Eq,-,  {8  ^ Short  of  this,,  which  is  an 
impractical  task  for  most  practical  hardware  sy  items,  a maximum  Iikli 

--W 

hood  receiver  can  b<  designed  to  recognize  the  peak  in  the  W(T„k  «t  |i 
curve  plotted  ae  a fcnttion  of  7 . k and  tq  . if  the  curve  is  symmetric 
about  its  peak  that  r lad  of  values  (7*  ko„  r.  corresponding  to  the  peal 
will  be  identical  to  the  triad  of  expectation  value  s (<  t > „ <kQ>  - <’T0>  and 
the  raaodrauir-  likli'hood  receiver  will  provide  the?e  expectation  values  a » its 

output  A further  important  assumption  is  usually  made  that  the  a prici  i 
prooab?lxiy  iV  it  k i i is  siowly  varying,  as  compared  with  the  wore 
sharply  p:!ak?.d  uf  '7,k  »t  ) „ or  Wff  jT,  k , t ) , This  assumption  is 
valid  when  time  base  distortion!  and  hence  the  r fa*  ter,  do  not.  exist 
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bn  a the  present  cas-  it  la  exactly  the  :j  • tht:  w 

upon  the  location  of  he  <»|1  o'  V>f  t it  : *’h,«  h cannot  be  ignored 

that  lei  ’ " o '.r  „ , iuriim  • hat  anpptinsum  ~f  e>  /«»r  is  ar  adaptive  re 

a v*;r  nihei  tha  - ' ordinary  coherent  or  -el-  t > i «*  stiAi  true  that 

iV  c . t , th,*  * l.  ic-ffj  probability  on  5ocat.oa  anc  ired  of  the  target  is 
o o o 

slowly  vary:, . g o'  er  that  range  oi  h ar d t to  * t oh  the  remainder  of 
the  right  fu.de  of  Kq  • 5 var  e.  r gnif  n antly  a-u  1 « 'em:  o ill  be  ignored 
ir>  establishing  the  peak  of  Eq  8 

■V e oust..  1 i -•  i 1 the  oe  ik  of 


W |t1  Ltf  i*.l  , r or 
o o o 


'V'  W(T|T,k.t 
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where  we  have  used  Eq.  17’  for  VV  (t  si  c Eq  for  W'f  It  , k t i , 

11  o 1 o o 

since  W(f It  ,2t  , r W {»  >y  virtue  of  Eq.  U 

* o o c 

i j v ■ o'.tsl y I-q  (9*  a •.  ains  pet  k vain  for  tn*  nnunt  positive  lvalue 
o;'  the  erporent.  Revritng  <9i. 


The  tight  side  of  Eq»  (11)  would  be  the  ordinary  correlation  integral  if  all 

were  set  equal  to  zero  and  the  second  term  were  missing,  corresponding 

i 

to  a complete  lack  of  time  base  distortions,,  In  the  presence  of  time  base 
distortions?  it  is  necessary  to  maximize  the  sura  of  the  two  terms  on  the 
right  side  of  Eq.  (11)  to  define  an  optimum  receiver  Obviously  the  first 
term,  or  correlation  integral?  should  be  made  as  large  as  possible  con- 
sistent with  a small  value  for  the  second  term?  A normal  procedure  in 
this  direction  would  be  to  arrive  at  a proper  choice  for  and  k . the 
average  delay  and  time  base  distortion  factors?  and  then  to  optimize  the 
correlation  integral  by  adjusting  the  instantaneous  distortion  factors  *t 
while  at  the  same  time  keeping  the  second  term  suitably  small,,  If  the  sec 
ond  term  were  not  present,,  it  would  be  possible  to  match  exactly  the  a.  to 
the  f . in  the  correlation  integral  by  adjusting  the  r.  ? but  reference  to 

i a 

Eq.  (1  would  show  that  this  procedure  merely  adapts  the  signal  components 

8,  to  the  noise  ? which  also  has  2WT  independent  components.  However. 

the  presence  of  the  second  term  in  Eq„  (11)  prevents  physically  unlikely  ad 

justments  of  the  t\  , In  particular?  this  term  will  become  large,  giving  a 

small  value  to  W <r  ,,  Eq„  |7i?  whenever  the  r are  varied  so  rapidly  that 
o i 

the  Tit1*  bandwidth  W is  exceeded  Thus?  the  correlation  matrix  (p  ? 

T 

which  in  the  case  of  sonar  is  determined  by  sea  conditions  and  targe',  accelera 
tion  capability,  will  set  restraints  upon  the  choice  of  a 2WT  set  of  components 
r.  to  maximize  Eq,  (III.  In  particular?  the  rapidity  with  which  the  t.  are 
varied  will  be  established  by  the  correlation  matrix  (p!  ? and  the  amplitude 

of  the  t variations  will  be  indicated  by  the  variance  <r  J"  , which  ia  also 

l 

a function  of  sea  state  and  target  motion.  When  the  average  noise  power 

<r  ^ is  small,,  the  coefficient  of  th*  first  term  in  Eq,  (11)  is  large?  and  this 
n 

term  becomes  relatively  more  important  than  the  second  In  this  high  signal 
to-noAse  case  it  is?  therefore?  legitimate  to  allow  the  to  assume  oven 
relatively  unlikely  sequences  of  values  in  the  interest  of  maximizing  the 


first  term,  even  though  the  second  term  becomes  more  negative  This 
merely  reflects  the  fact  that  when  the  noise  power  is  low?  there  is  a rela 
tively  high  degree  o{  confidence  that  whatever  time  base  distortions  occur 
in  the  signal  are  produced  by  the  medium  and  not  the  additive  noise,  despite 
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the  lower  a priori  probability  which  may  be  associated  with  this  type  of 
medium  distortion  Of  course,,  the  reverse  is  true  when  cr ^ is  small. 

indicating  relatively  calm  sea*  or  the  fact  that  targets  which  may  be  la 
the  surveillance  area  have  relatively  small  acceleration  capability  in  the 
signal  duration  T,  or  when  <r  n is  large,  indicating  the  presence  of  a large  noi 

contribution  In  this  case,  the  Second  term  w Eq„  (11)  becomes  more 

important  and  the  tolerable  departures  of  the  t.  from  their  zero  mean 

are  small.  These  considerations  lead  quite  ^.turally  to  certain  design 

forms  for  an  adaptive  receiver,  which  are  discussed  in  Technical  Memo- 
7 

randum  No.  65. 


1 Ross  E Williams.  ‘"The  Design  of  Adaptive  Receivers  for  Optimum  De 
tection  of  Temporally  Distorted  Signals"  (Hudson  Laboratories,  Columbia 
University  Technical  Memorandum  No.  65,  January  14,  1963). 
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